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Abstract 

Inclusive jet differential cross sections for the reaction e + p — ► e + + jet + X with quasi- 
real photons have been measured with the ZEUS detector at HERA. These cross sections 
are given for the photon-proton centre-of-mass energy interval 134 < W < 277 GeV and 
jet pseudorapidity in the range — 1 < rf et < 2 in the laboratory frame. The results are 
presented for three cone radii in the r/ — <p plane, R = 1.0, 0.7 and 0.5. Measurements of 
da/dr] jet above various jet-transverse-energy thresholds up to 25 GeV and in three ranges 
of W are presented and compared to next-to-leading order (NLO) QCD calculations. For 
jets defined with R = 1.0 differences between data and NLO calculations are seen at high 
rf et and low E^ . The measured cross sections for jets defined with R = 0.7 are well 
described by the calculations in the entire measured range of rf et and E^. The inclusive 
jet cross section for E 3 ^ 1 > 21 GeV is consistent with an approximately linear variation 
with the cone radius R in the range between 0.5 and 1.0, and with NLO calculations. 
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1 Introduction 



At HERA, photon-proton reactions are studied by means of ep scattering at low four-momentum 
transfers squared (Q 2 ~ 0). In photoproduction, two types of QCD processes contribute to 
the production of jets |],|2|] at leading order (LO): either the photon interacts directly with 
a parton in the proton (the direct process) or the photon acts as a source of partons which 
interact with those in the proton (the resolved process). Differential cross sections for inclusive 
jet photoproduction using a cone algorithm have been previously presented as a function of 
the jet pseudorapidityj] (r/ je ') and transverse energy (E^) for Ej? up to 17 GeV ||f|,[J. The 
calculated cross sections depend on the proton parton distributions in the region of Bjorken-x 
above approximately 10~ 2 , where they are well constrained by other measurements ||. Such 
jet measurements therefore offer a potential means of studying the parton distributions in the 
photon [M0,p,^[rn,|nj at higher scales than those probed in e + e~ interactions |T3|. However, 
various aspects of the comparison between theory and experiment need to be addressed before 
a reliable determination of the photon parton distributions can be made. 

Next-to-leading order (NLO) calculations including resolved plus direct processes and using 
the NLO parametrisations of the photon parton distributions of GRV [13 have been compared 
H to our previous measurements M . Discrepancies were observed in the forward region {rf et > 
1) for low E^ 1 [E^ 1 ~ 8 GeV) which prevented any strong conclusion being drawn on the 
photon parton distributions. Moreover, the comparison of the transverse energy flow between 
data and leading-logarithm parton-shower Monte Carlo simulations [|J showed a discrepancy 
at high rf et that could be attributed to energy not associated with the hard-scattering process 
(the 'underlying event'). Such an underlying event is not included in the NLO calculations 
and, therefore, the comparison between data and NLO calculations becomes problematic. The 
transverse energy inside the cone of the jet in the r\ — (p plane due to the underlying event is 
naively expected to be proportional to the area covered by the cone. Therefore, measurements 
performed with different cone radii can elucidate the effects of a possible underlying event. In 
addition, the NLO calculations for jets defined with a cone radius R ~ 0.7 are expected to 
be most stable with respect to variations of the renormalisation and factorisation scales ||14|| . 
Measurements of the jet cross sections in different ranges of the 7p centre-of-mass energy (W) 
provide a further means of comparing data and calculations. 

In this paper, measurements of da/drf et are presented for various cone radii. In each case, 
measurements of da/drf et integrated above four different E^ thresholds (14, 17, 21 and 25 
GeV) are performed. For R = 1.0, the measurement of da/drf et is extended to higher E^ 1 
values as compared to the previous data ||. We have, in addition, a better understanding of 
the energy scale of the jets. First measurements of do jdrf et in three regions of W are presented 
for R = 1.0 and 0.7. The dependence of the inclusive jet cross section on the jet cone radius is 
presented. NLO calculations P,|IP| which include resolved plus direct processes are compared 
to the measurements. 

The data sample used in this analysis was collected with the ZEUS detector in e + p interac- 
tions at the HERA collider and corresponds to an integrated luminosity of 2.65 pb _1 , which is 
a five- fold increase in statistics over the previous analysis pj. 



1 The ZEUS coordinate system is defined as right-handed with the Z axis pointing in the proton beam 
direction, hereafter referred to as forward, and the X axis horizontal, pointing towards the centre of HERA. 
The pseudorapidity is defined as r\ = — ln(tan|), where the polar angle 9 is taken with respect to the proton 
beam direction. 
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2 Experimental conditions 



During 1994 HERA operated with protons of energy E p = 820 GeV and positrons of energy 
E e = 27.5 GeV. The ZEUS detector is described in detail in 116]] - The main subdetectors 
used in the present analysis are the central tracking system positioned in a 1.43 T solenoidal 
magnetic field and the uranium-scintillator sampling calorimeter (CAL). The tracking system 
was used to establish an interaction vertex and to cross-check the energy scale of the CAL. 
Energy deposits in the CAL were used in the jet finding and to measure jet energies. The CAL 
is hermetic and consists of 5918 cells each read out by two photomultiplier tubes. Under test 
beam conditions, the CAL has energy resolutions of 18%/^/E for electrons and 35%/\/E~ for 
hadrons. Jet energies are corrected for the energy lost in inactive material in front of the CAL 
which is typically about one radiation length (see Section 5). The effects of uranium noise 
were minimised by discarding cells in the inner (electromagnetic) or outer (hadronic) sections 
if they had energy deposits of less than 60 MeV or 110 MeV, respectively. The luminosity was 
measured from the rate of the bremsstrahlung process e + p — *> e + p^. A three-level trigger was 



used to select events online [16, 17] 



3 Data selection and jet search 

Offline, events from quasi-real photon-proton collisions were selected using similar criteria as 
employed previously Q. The main steps are briefly discussed here. The contamination from 
beam-gas interactions, cosmic showers and beam-halo muons is negligible after demanding: 
a) at least two tracks pointing to the vertex; b) the vertex position along the beam axis to 
lie in the range — 29 < Z < 36 cm; c) fewer than five tracks not associated with the vertex 
and compatible with an interaction upstream in the direction of the proton beam; and d) the 
number of tracks not associated to the vertex be less than 10% of the total number of tracks. 
Deep-inelastic (DIS) charged-current e + p scattering events are rejected by requiring the total 
missing transverse momentum (p/?) to be small compared to the total transverse energy (Etf*)'- 
j/t I v/ 1 Eip 1 < 2 GeV^ . DIS neutral-current events with an identified scattered positron candidate 



in the CAL, according to the algorithm described in ||18|| , are removed from the sample. The 
selected sample consists of events from e + p interactions with Q 2 < 4 GeV 2 and a median of 
Q 2 fa 10~ 3 GeV 2 . The events are restricted to the kinematic range 134 < W < 277 GeV using 
the procedure described in Section 5. 

An iterative cone algorithm in the rj — (p plane [0, |2(| is used to reconstruct jets from 
the energy measured in the CAL cells. A detailed description of the algorithm can be found 
in [|17]]. The jets reconstructed from the CAL cell energies are called cal jets and the variables 
associated with them are denoted by E^lai, rfj^ and ip 3 ^. The axis of the jet is defined according 



to the Snowmass convention ||20|| , where rf^ (v^) is the transverse-energy weighted mean 
pseudorapidity (azimuth) of all the CAL cells belonging to that jet. Events with at least one 
jet satisfying E^ 1 ^ > 10 GeV and — 1 < rf c ^ < 2 are retained. Three samples of jets have 
been selected depending on the cone radius used in the jet search: 18897 jets for R = 1.0, 
11197 jets for R = 0.7 and 7070 jets for R = 0.5. The only significant remaining background is 
from unidentified DIS neutral current interactions with Q 2 > 4 GeV 2 , which is estimated using 
Monte Carlo techniques to be below 2%. 



2 



4 Monte Carlo simulation 



Samples of events were generated using Monte Carlo (MC) simulations to determine the re- 
sponse of the detector to jets of hadrons and the correction factors for the inclusive jet cross 
sections. 



The programs PYTHIA 5.7 JzlJ] and HERWIG 5.8 |22j were used to generate photoproduction 
events for resolved and direct processes. In PYTHIA the positron-photon vertex was modelled 
according to the Weizsacker- Williams approximation. In the case of HERWIG, the exact matrix 
elements were used for direct processes (e + g — > e + qq and e + q — > e + qg) and the equivalent 
photon approximation for resolved processes. Events were generated using GRV-HO [13| for the 



photon parton distributions and MRSA [23| for the proton parton distributions. In addition, 
samples of events using the LAC1 parametrisation [24] for the photon parton distributions 
were considered. In both generators, the partonic processes were simulated using LO matrix 
elements, with the inclusion of initial and final state parton showers. Fragmentation into 
hadrons was performed using the LUND [2~5| string model as implemented in JETSET [BBJ in 
the case of PYTHIA, and the cluster model [27j in the case of HERWIG. Samples of events were 
generated with different values of the cut-off on the transverse momentum of the two outgoing 
partons starting at pTmin = 8 GeV. For the measurements presented in this paper, the events 
generated using the PYTHIA and HERWIG programs have been used for calculating energy 
corrections and for correcting the data for detector and acceptance effects. The corrections 
provided by the PYTHIA generator have been used as default values and the ones given by 
the HERWIG generator have been used to estimate the systematic errors coming from the 
fragmentation model. 

Additional samples of events were generated using the option of multiparton interactions 
(MI) in PYTHIA. This option, which applies only to resolved processes, adds interactions 
between the partons in the proton and the photon remnants to the hard scattering process of 
the event. These multiparton interactions are calculated as LO QCD processes and give an 
estimation of the underlying event. The PYTHIA MI events were generated with a cut-off 



for the effective minimum transverse momentum for multiparton interactions of 1 GeV [1/ 



and with a cut-off on the transverse momentum of the two outgoing partons from the hard 
scattering of PTmin = 8 GeV. 

All generated events were passed through the ZEUS detector and trigger simulation programs 
rB| . They were reconstructed and analysed by the same program chain as the data. 



For the Monte Carlo events, the jet search is performed from the energy measured in the CAL 
cells in the same way as in the data. The same jet algorithm is also applied to the final state 
particles. In this search, all particles with lifetimes longer than 10~ 13 s and with polar angles 
between 5° and 175° are considered. The jets found are called hadron jets and the variables 
associated with them are denoted by -E^Ld; f]j^ d , and fj^d- Hadron jets with E^ had > 14 GeV 
and — 1 < rf^ d < 2 are selected. 



5 Energy corrections 



The fivefold increase in statistics in 1994 allowed the CAL energy scale to be studied in more 
detail than in ref. 0]. The comparison of the energy measured in the central region of the CAL 
to the momentum measured in the tracking system for the scattered positron in neutral current 
DIS events, and the transverse momentum balance in neutral current DIS events, showed a 



(6 ± 3)% difference between data and MC [P8| . This 6% disagreement has been corrected for 
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in the present analysis. In the analysis of the 1993 data, the possibility of such a discrepancy 
was allowed for in the systematic uncertainties. 

Particles impinging on the CAL lose energy in the inactive material in front of the CAL. 
The inactive material constitutes about one radiation length except in the region around the 
rear beampipe, 9 > 170°, and the support structures, 25° < 9 < 45° and 130° < 9 < 145°, where 
it reaches 2.5 radiation lengths. For the measurements presented here, the transverse energy of 
the jets has also been corrected for these energy losses as explained below. 

The comparison of the reconstructed jet variables between the hadron and the cal jets in 
simulated events shows no significant systematic shift in the angular variables if^ and with 
respect to rf^ and tpj^d- Therefore, no correction is needed for "if et and ip jet {if et ss 77^* and 
ipi et pa v?ca/)- However, the transverse energy of the cal jet underestimates that of the hadron jet 
by an average amount of 16% with an r.m.s. of 11%. The transverse energy corrections to cal 
jets averaged over the azimuthal angle were determined using the MC events. These corrections 
are constructed as multiplicative factors, C(E^ cal ,r] 3 ca \), which, when applied to the Et of the 
cal jets provide the 'corrected' transverse energies of the jets, E^ 1 = C(E^ t ccd ,rf ( ^ l ) x E^lai- 

The method of Jacquet-Blondel |2"9|| , applied to the photoproduction regime [BO], is used to 



estimate W from the energies measured in the CAL cells: W = y2E p ■ {E — pz), where E is 
the total CAL energy and pz is the Z component of the directed energy measured in the CAL 
cells. Due to energy lost in the inactive material in front of the CAL and to particles lost in the 
rear beampipe, W cal systematically underestimates W by approximately 10% with an r.m.s. of 
5%. This effect is adequately reproduced by the MC simulation of the detector. To compensate 
for this underestimation, MC samples of events were used to determine a correction procedure 
to W cal as a function of W cal and of the pseudorapidity of the most backward jet in the event 
iVmin)- This correction has been constructed as a multiplicative function, Y(W cal , rj^ n ), in 
a similar way as the correction to the jet transverse energy. When applying the function Y 
to W cal , W = Y(W cal , r/min) x W cal , the corrected 7p centre-of-mass energy is obtained, and 
events with 134 < W < 277 GeV are retained. 



The response of the CAL to jets has been checked by the following procedure pl| . In 
the central region (jr/- 76 *! < 1), the multiplicity distribution and the p^-spectrum of charged 
particles within the cal jets have been compared for data and Monte Carlo samples using 
the reconstructed tracks. The tracks were required to be in the ranges \rj track \ < 1.5 and 
ptrack > 3QQ MeV, where p t f ack is the transverse momentum of the track with respect to the beam 
axis. Tracks were associated with a cal jet when the extrapolated trajectory reached the CAL 
within the cone of the cal jet. PYTHIA describes well all the measured distributions. In this 
rf et region, the momenta of the tracks in the cal jet are used to determine the total transverse 
energy carried by the charged particles, EjT tracks . Then, the ratio r tracks = E^ racks /E^ cal is 
formed, and the distributions of this ratio for the inclusive cal jet sample with R = 1.0 in 
data and simulations are compared, as shown in Figure [TJa. The mean value of the distribution 
in r t racks has been determined as a function of rf et for data (<r tra cks>data) and simulations 
tracks^" mc)- From the values of the quantity {<r tracks > data / <r tr acks>Mc) ~ L shown in 
Figure [TJc (circles), we conclude that the energy scale of the jets with |7y je '| < 1 is correct to 
within ±3%. 

In the forward region, 1 < n^ et < 2, the energy scale of the jets is studied using the transverse 
energy imbalance in dijet events with one jet in the central region and the other in the forward 
region. The distributions of the ratio r di j et = E^l al (forward jet) /E^l al (central jet) in data and 
simulations are compared in Figure |TJh>. The values of the quantity (<rdij e t>data / < r dijet>Mc) — l 
(see Figure 0c, square symbols) show that in the forward region the energy scale of the jets is 
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also correct to within ±3%. 

It is noted that since the widths of the r tracks and r di j et distributions in the data are reasonably 
well described by the PYTHIA simulations, the resolution in the energy of the jets is also 
correctly described. 

This procedure has been also applied to the inclusive cal jet sample with R = 0.7 and leads to 
the same conclusions. The use of HERWIG instead of PYTHIA gives similar results. Therefore, 
a ±3% uncertainty on the energy scale of the jets is included as a systematic variation in the 
present analysis. 



6 Jet profiles 

The presence of energy not associated to the hard-scattering process (the 'underlying event') 
in the data has been investigated through the study of the transverse energy flow around the 
jet axis both inside and outside of the jet cone. 

The transverse energy profile around the jet axis was measured using the energies and angles 
of the CAL cells uncorrected for detector effects. The distribution of transverse energy in the 
hemisphere of the jet, as a function of Arj = r\ ce ii~rf et and integrated over |A<^| = \(p ce u — f 3et \ < 
7r/2, is shown in Figures |2] and [3] for the inclusive jet data samples (R = 1.0 and 0.7) in three 
^jet ran g es anc i t wo E 3 ^ 1 regions^ The data exhibit a pronounced peak at Arj = and an 
asymmetric pedestal. The height of the peak increases as E 3 ^ 1 increases. As a function of 
r) jet , it is fairly constant in the region — 1 < rf et < 1 and decreases in the region r] jet > 1; this 
decrease is most significant for jets with R = 1.0. The height of the pedestal for At] > 1 (proton 
side) slightly increases with increasing r] jet , the effect being more pronounced for R = 1.0 and 
low Ej? . 

The expectations from PYTHIA simulations including resolved plus direct processes are 
compared to the data in Figures ^| and The transverse energy profile in the data is well 
described by the simulations of PYTHIA except for jets with rj jet > 1 and lowest Ej? (E 3 ^ 1 m 
14 GeV). In this region, an excess of transverse energy outside of the jet cone with respect to 
PYTHIA simulations is observed P,^|,p|,p2|. The excess is reduced for jets defined with R = 0.7 
in comparison to jets defined with R — 1.0 (see Figure §). In order to simulate an increased 
energy flow, the PYTHIA MI generator is used, which gives rise to energy not associated with 
the hard-scattering process. PYTHIA MI gives an improved description of the data for forward 
low-E^ jets with R = 1.0, but lies above the data for r] jet < 1 in the case of R = 1.0 and 
in all r) 3et ranges for R = 0.7. For jets with Wf > 21 GeV, no significant discrepancies are 
observed between data and PYTHIA simulations with or without multiparton interactions (see 
Figure |). 

The internal structure of the jets may be investigated using the jet shape, defined as the 
average fraction of the jet transverse energy that occurs inside an inner cone concentric with 
the jet defining cone fll4]| . The shape of jets selected using R = 1.0 has been recently measured 



in photoproduction at HERA [H]] and found to be well described by the PYTHIA calculations 
except for the inclusive production of jets with if et > 1 and low E 3 ^ 1 (14 GeV < E 3 ^ < 
17 GeV). We have performed the same type of analysis for jets with R = 0.7 and, in this 
case, the measured jet shapes (not shown) are well described by the PYTHIA (with or without 
multiparton interactions) calculations in the entire if et region. 

2 The decrease of the A77 distribution seen both in data and the simulations in the region A?7 > 2 for the 
forward jets (1 < T\> et < 2) is a geometric effect: the most forward edge of the CAL is at r\ — 4.3. 



5 



These observations indicate that the uncertainties on the jet measurements due to possible 
underlying event contributions become reduced at high E 3 ^ 1 [E 3 ^ 1 > 21 GeV) or when using a 
reduced cone radius (R = 0.7). 

7 Acceptance corrections and systematic uncertainties 

The MC generated event samples of resolved and direct processes were used to compute the 
acceptance corrections to the inclusive jet distributions. These correction factors take into 
account the efficiency of the trigger, the selection criteria and the purity and efficiency of the 
jet reconstruction. The differential cross sections da/di] jet are then obtained by applying bin- 
by-bin corrections to the measured jet distributions. The predictions of the generators PYTH1A 
and HERWIG for the uncorrected distributions were compared to the data for several choices 
of the parton densities in the photon and proton and for various combinations of resolved and 
direct processes. A good description of the rf et data distributions is obtained by the MC except 
for forward low-E^ jets with R = 1.0. The bin-by-bin correction factors lie between 0.7 and 
1.4 depending on rj jet , E 3 T et threshold and W region considered. The dominant effect arises from 
migrations over the E^ et threshold. 

A detailed study of the sources contributing to the systematic uncertainties of the measure- 
ments has been performed. The study of the systematic uncertainties includes (a typical value 
for each item is indicated): 

• Use of the HERWIG generator to evaluate the energy corrections to cal jets and the 
correction factors to the observed inclusive jet distributions. The effect of this variation 
is typically within ±5% in the region 0.5 < rf et < 2 and increases to 10% for r] jet < 0.5. 

• Uncertainties in the simulation of the trigger and the variation of the cuts used to select 
the data within the ranges allowed by the comparison between data and MC simulations 
(« 5%). 

• Use of the PYTHIA generator including multiparton interactions in resolved processes 
to evaluate the energy corrections to cal jets and the correction factors to the observed 
inclusive jet distributions (~ 3%). In the region of forward low-E^ 1 jets with R = 1.0, 
an improved description of the data is obtained by using PYTHIA MI. 

• Choice of different parton densities in the photon (GRV-HO and LAC1) for the generation 
of the PYTHIA MC samples (« 2%). 

All these systematic uncertainties have been added in quadrature to the statistical errors and 
are shown as thin error bars in the figures. 

• The absolute energy scale of the cal jets in simulated events has been varied by ±3% 
for the reasons discussed in Section 5. The effect of this variation on the inclusive jet 
cross sections is ~ ±12% in the region < i] jet < 2, and increases up to ~ 35% for 
rf et « —1. This uncertainty represents the dominant source of systematic error and is 
highly correlated between measurements at different rf et points. It is shown as a shaded 
band in each figure. 

In addition, there is an overall normalisation uncertainty of 1.5% from the luminosity deter- 
mination, which is not included. 
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8 Results 



8.1 Differential cross sections 

We present measurements of inclusive differential jet cross sections for the reaction 

e +p _> e + 4. j e t + x 

in the kinematic region defined by Q 2 < 4 GeV 2 and 134 < W < 277 GeV. These cross sections 
refer to jets at the hadron level with cone radii of R = 1.0 and 0.7 units in the 77 — (p plane. The 
cross section da/drf et has been measured in the rf et range between —1 and 2 integrated above 
£j. e< from four different thresholds {E 3 ^ > 14, 17, 21 and 25 GeV). The cross section da/drf** 
for Eq? > 14 GeV has also been measured for three different regions of W: 134 < W < 190 
GeV, 190 < W < 233 GeV and 233 < W < 277 GeV. The results are presented in Figures | 
to and in Tables [l] to 

For E^* > 14 and 17 GeV, the behaviour of the cross section as a function of rf et in the 
region rf et > 1 is very different for R — 1.0 and R — 0.7 (see Figures |4] and ||): it is constant for 
R = 1.0 whereas it decreases as rf et increases for R = 0.7. On the other hand, the behaviour 
for E 3 ^ 1 > 21 and 25 GeV is approximately the same in both R = 1.0 and 0.7 cases. There 
are two effects which contribute to the observed differences in the region rj ]et > 1 for E 3 ^ 1 > 14 
and 17 GeV: a) the jets become broader as r] jet (E^) increases (decreases) [|7]] and, b) the 
pedestal in the jet profile is integrated over approximately half the area for jets with R = 0.7. 
In addition, the height of the pedestal (see Section 6) is larger for forward jets with R = 1.0. 
Therefore, the differences between the cross sections for the two radii can be attributed to the 
fact that the use of R = 0.7 selects more collimated jets and suppresses the underlying event 
contribution. 

The results for da/df] jet in different regions of W for E 3 T et > 14 GeV and with R = 1.0 
(R = 0.7) are presented in Figure ^ (0). For R = 1.0, the cross section increases with increasing 
values of if et and is constant in the high if et region, whereas for R = 0.7 the cross section 
decreases as if et increases in the high t] jet region. For increasing values of W the maximum of 
the cross section with R = 1.0 (R = 0.7) shifts to lower values of i] jet . As the energy of the 
incoming quasi-real photon increases, W increases and the events are boosted more backwards 
in the laboratory frame. 



8.2 Comparison to NLO calculations 



NLO QCD calculations of dajdrf ei are compared to our measurements in Figures § to [7]. 

These predictions include resolved and direct processes. The CTEQ4M |33j proton parton 
densities have been used. For the photon parton distributions, the AFG |34| , GRV-HO |1^] and 
GS96 PS| parametrisations have been usedQ. In the calculations shown here, the renormalisation 
and factorisation scales have been chosen equal to E^ 1 and a s was calculated at two loops with 
A^ = 296 MeVi. 

The comparison of the data with NLO calculations is subject to the uncertainty in matching 
the experimental and theoretical jet algorithms. Since the calculations include only up to three 
partons in the final state, the maximum number of partons in a single jet is two. Therefore, the 
overlapping and merging effects of the experimental jet algorithm are not reproduced in the 



The calculations using GRV-HO or GS96 are from [Hi and those using AFG from (TJ. For the same photon 
parton distributions, the calculations from || and [jlfj differ typically by less than ±5%. 
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theoretical calculation [T^, 36]. An attempt was made to simulate these effects by introducing 



an ad-hoc Rsep parameter fl4fl : two partons are not merged into a single jet if their separation 
in the r] — ip plane is more than Rsep- The calculations of the cross sections shown in Figures f| 
to |7| have been made for Rsep = R- In addition, the calculations using GS96 and Rsep = 2-R 
are also shown. The spread of the calculations using GS96 for Rsep = R and Rsep = 2R 
indicates the magnitude of the theoretical uncertainty due to these effects. 

As discussed above, the NLO calculations refer to jets built out of at most two partons 
whereas the measurements refer to jets at the hadron level. An estimate of the effects of hadro- 
nisation has been obtained by comparing the cross sections for jets of hadrons and jets of partons 
calculated with the PYTHIA generator. The ratio of (G?cr/(i?7 : ' e *[hadrons])/((io"/(i?7- 7e '[partons]) 
for jets with R = 1.0 (R = 0.7) is relatively constant as a function of r] jet and within ap- 
proximately 10% (20%) of unity. Due to the approximations used in the MC simulations, 
these estimations are not to be taken as corrections to the parton level for the measurements 
presented here. 

The NLO calculations give a good description of the measured differential cross sections in 
magnitude and shape for E^* > 21 and 25 GeV for both cone radii R = 1.0 and 0.7. For 
Ej? 1 > 14 GeV, the behaviour of the measured cross sections is different for R — 1.0 and 0.7, 
whereas the calculations exhibit the same shape for both radii. For R = 1.0, the shape of the 
cross section is well described for —1 < rf et < 0.5. For higher values of r] jet , the measured 
cross section is constant, as discussed in Section 8.1, whereas the theoretical curves decrease. 
These differences are not present when R = 0.7 is used: the NLO calculations describe well 
the magnitude and shape of the measured differential cross sections with R = 0.7 for all E 3 ^ 1 
thresholds in the entire range of rf et . 

The measured differential cross sections in ranges of W for R = 1.0 and E 3 ^ 1 > 14 GeV are 
reasonably well described for low values of r/ je *, whereas the NLO calculations fail to describe 
the high rf et region. The excess of the measured cross section with respect to the calculations 
increases with increasing W. On the other hand, the measured differential cross sections in 
bins of W are reasonably well described by the NLO calculations using R = 0.7 in the entire 
region of rf et . 

The failure of the NLO calculations to describe the measured cross section for forward low- 
E J T et jets with R = 1.0 may be due to the following effects: a) the uncertainty due to the choice 
of renormalisation and factorisation scales is larger than in the case R = 0.7 (see next section), 
and b) non-perturbative contributions like that of the underlying event, which is reduced for 
jets with R = 0.7, are not included. On the other hand, for jets defined with R = 0.7 the 
measured cross sections are well described by the calculations and the uncertainties on the 
measurements are comparable to the spread of the predictions using different parametrisations 
of the photon parton distributions. 



8.3 Cone radius dependence of the cross section 

The cone radius dependence of the inclusive jet cross section, a(R), has been studied. Measure- 
ments have been performed of the inclusive jet cross section integrated above E 3 ^ 1 > 21 GeV 
and —0.5 < r] jet < 2 for three different cone radii (R = 1.0, 0.7 and 0.5). These cross sections 
are given in the same Q 2 and W kinematic region as the measurements presented in Section 8.1. 
As observed in the jet profiles (see Section 6), the uncertainties on the jet cross sections due to a 
possible underlying event become reduced at E^ 1 > 21 GeV. The results for a(R) are presented 
in Figure || and Table ||. The measured cross section is consistent with a linear variation with 
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R in the range between 0.5 and 1.0. 

The results of LO and NLO QCD calculations of cr(R) @, which are performed at the parton 
level, for different values of the renormalisation and factorisation scales [i are shown in the inset 
of Figure |[ The LO and NLO GS96 (CTEQ4) sets of photon (proton) parton densities have 
been used. The LO predictions do not depend on R since there is only one parton per jet and 
show a large variation with /i. NLO calculations give the lowest-non-trivial order independent 
contributions to the jet cross section and the \i dependence is largely reduced. However, at 
small (large) values of R, the NLO predictions for cr(R) become a monotonically increasing 
(decreasing) function of fx. The calculations are most stable for R rs 0.5 — 0.7, consistent with 



the conclusions of JT4j. The uncertainty on the predicted cross section due to the choice of fx, 
estimated by changing ll from E 3 T et /4 to E J T Gt , is 5% (20%) at R = 0.7 (R =1.0). 

The slope of cr(R) depends on the choice of fx, and is largest (smallest) for small (large) 
values of fx (see inset of Figure |8|). The slope of c(R) in the NLO calculation with fx = E^/A 
is closest to that of the measured cross section. In addition to the uncertainty coming from 
the choice of fx, the predictions are affected by the value of Rsep- QCD calculations with 
ji = Ej? /A and for two values of Rsep, Rsep = R and 2R, are compared to the measurements 
in Figure §. Since the LO predictions of the inclusive jet cross section do not depend on R, the 
data show the need for QCD corrections. The NLO calculations are consistent with the data 
within the theoretical and experimental uncertainties, both of which are at the 20% level. 



9 Summary and conclusions 

Measurements of differential cross sections for inclusive jet photoproduction in e + p collisions at 
HERA using the data collected by ZEUS have been presented. The cross sections refer to jets 
at the hadron level found with an iterative cone algorithm in the f] — <f plane. Measurements 
of the jet cross sections with two different cone radii, R = 1.0 and 0.7, have been performed. 
These cross sections are given in the kinematic region defined by Q 2 < 4 GeV 2 and 134 < W < 
277 GeV. 

A comparison has been made of the transverse energy profiles around the jet axis between 
data and the leading-logarithm parton-shower simulations of PYTHIA. Requiring high Ej? 
(Ej? > 21 GeV) or using a cone radius of R = 0.7 reduces the discrepancy between data and 
PYTHIA in the forward region. 

NLO QCD calculations |||TIJ using currently available parametrisations of the photon parton 
distributions are compared to the measured cross sections. The uncertainties on the calculations 
due to the choice of renormalisation and factorisation scales, and non-perturbative effects like 
the underlying event are smaller for jets with R = 0.7 than in the case of R = 1.0. The 
calculations describe the measured cross sections well for jets defined with R = 1.0 and 0.7 for 
E^ > 21 and 25 GeV. At lower values of E^ 1 differences between data and the calculations are 
seen in the forward region for jets defined with R = 1.0. On the other hand, the calculations 
describe well the measured differential cross sections in the entire range of rf et for jets defined 
with R = 0.7. These conclusions are reinforced when the data are considered in different ranges 
of W. The uncertainties on the measurements with R — 0.7 are comparable to the spread of 
the predictions using different parametrisations of the photon parton distributions. 

The measured cross section for jets with E^ 1 > 21 GeV and —0.5 < rf et < 2 is consistent 
with a linear variation with the cone radius R in the range between 0.5 and 1.0, and shows the 
need for QCD corrections. The NLO calculations are consistent with the data within the 20% 
theoretical and experimental uncertainties. 
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0.62 


1790 ± 50 ± 60 


(+220, —160) 


n 88 

U.oo 


1 ( OO It ou It oU 


f-Lonn 1 p.c\\ 
^-i-zuu, — lOU J 


l.lz 


1 ( 10 ± OU ± 1 1U 


(_ + ! ( U, — 14U J 


1 QS 
l.OO 


i «on -l K(\ -i- an 
loyU ± OU ± oU 


( i i on 1 7n"\ 

^ + loU, — 1 ( UJ 


1 fiO 
l.OZ 


1000 it OU it 11U 


^-|-zoU, — 10U J 


1 88 
l.OO 


I 78 K 4- KC\ 4- 1 nn 

I I oO it OU it 1UU 


f-LOon oin"\ 
^-hzzU, — ziu J 




> 17 GeV 




-0.62 


80 ± 10 ± 10 


(+30, -20) 


-0.38 


185 ± 15 ± 30 


(+50, -40) 


-0.12 


355 ± 25 ± 30 


(+70, -50) 


0.12 


500 ± 25 ± 30 


(+80, -70) 


0.38 


625 ± 30 ± 50 


(+100, -90) 


0.62 


wi» p» i or i A r\ 

755 ± 35 ± 40 


(+100, -70) 


U.oo 


i ou it oo it zu 


(4-i nn an"\ 
^-|-1UU, — oU J 


110 

l.lz 


70 k -U ik -U c;n 
l zO ± oO ± OU 


1^+oU, — yuj 


1 QS 
l.OO 


«on -l Qn 4- /in 
0»U ± oU ± 4U 


( i on 7n"\ 
l^+yu, — /uj 


1 fiO 
l.OZ 


71 n -i- ik _i_ 7n 

( 1U it oO it iU 


('-Li nn Qn"\ 
^-|-iuu, — yu ) 


1 88 
l.OO 


««K _|_ on _|_ Qfl 
000 it oU it oU 


('-Li 1 n 7n"\ 
l^+iiu, — ( UJ 




E J T > 21 GeV 




n oc; 
-U.zo 


KC\ 4- K 4- 1 K 
OU it it 10 


('-Li 1 r\\ 
(,+J-O, — 1UJ 


U.zo 


one; _i_ m j- 1 k 
ZUO it 1U it 10 


l^+oo, — ooj 


n 7K 

U. l 


QC\K _|_ 1 K _|_ K 
OUO it 10 it 


('j-Cin /i 
l^+ou, — 4UJ 


1 oc; 

l.zO 


osn 4- i k 4- on 

ZoU It 10 It zu 


('-LQn ^^"i 
l^+ou, — OOJ 


1.75 


235 ± 15 ± 30 


f + 35 _35) 




E J T et > 25 GeV 




0.25 


70 ± 5 ± 5 


(+15,-10) 


0.75 


125 ±10 ±15 


(+20, -20) 


1.25 


125 ±10 ±15 


(+20,-15) 


1.75 


120 ±10 ±10 


(+25, -20) 



Table 1: Differential e + p cross section da/drf^ for inclusive jet production integrated above 
different E^ thresholds in the kinematic region defined by Q 2 < 4 GeV 2 and 134 < W < 
277 GeV for jets with a cone radius R = 1.0. The statistical and systematic uncertainties 
—not associated with the absolute energy scale of the jets— are also indicated. The systematic 
uncertainties associated to the absolute energy scale of the jets are quoted separately. The 
overall normalization uncertainty of 1.5% is not included. 
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rf et 


da/di]^ et ± stat. ± syst. [pb] 


syst. E^-scale [pb] 




E J T et > 14 GeV 




-0.88 


85 ±15 ±10 


(±30, -20) 


-0.62 


260 ± 20 ± 60 


(±70, -50) 


-0.38 


495 ± 30 ± 60 


(±110,-80) 


-0.12 


735 ± 35 ± 80 


(±120,-90) 


0.12 


935 ± 35 ± 50 


(±120,-110) 


0.38 


1050 ± 40 ± 60 


(±130,-120) 


0.62 


1225 ± 40 ± 60 


(±140,-120) 


0.88 


1165 ±40 ±40 


(±110,-110) 


1.12 


1065 ± 45 ± 70 


(±120,-110) 


1.38 


1055 ± 40 ± 30 


(±120,-100) 


1.62 


965 ± 35 ± 40 


(±110,-100) 


1.88 


850 ± 35 ± 40 


(±90, -90) 




E J T et > 17 GeV 




-0.62 


60 ±10 ±10 


(±20, -20) 


-0.38 


150 ±15 ±20 


(±30, -30) 


-0.12 


240 ± 20 ± 30 


(±40, -40) 


0.12 


380 ± 25 ± 40 


(±60,-50) 


0.38 


470 ± 25 ± 10 


(±80, -70) 


0.62 


525 ± 25 ± 20 


(±70, -70) 


0.88 


560 ± 30 ± 40 


(±70, -70) 


1.12 


495 ± 30 ± 30 


(±50, -50) 


1.38 


465 ± 25 ± 30 


(±50, -50) 


1.62 


395 ± 25 ± 40 


(±50, -40) 


1.88 


385 ± 25 ± 30 


(±40, -50) 




E J T et > 21 GeV 




-0.25 


35±5± 10 


(±10,-10) 


0.25 


160 ± 10 ±5 


(±30,-25) 


0.75 


210 ± 10 ±5 


(±30,-25) 


1.25 


195 ±10 ±10 


(±25,-25) 


1.75 


170 ± 10 ± 20 


(±25,-25) 




E 3 T et > 25 GeV 




0.25 


64 ± 7 ± 6 


(±10,-10) 


0.75 


93±8± 11 


(±15,-10) 


1.25 


85±8± 15 


(±15,-10) 


1.75 


85±8±8 


(±15,-10) 



Table 2: Differential e + p cross section dajdif et for inclusive jet production integrated above 
different E^ thresholds in the kinematic region defined by Q 2 < 4 GeV 2 and 134 < W < 277 
GeV for jets with a cone radius R = 0.7. Other details as in Table [l]. 
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da/dr]^ et ± stat. ± syst. [pb] 


syst. _E^ e *-scale [pb] 




134 < W < 190 GeV 




0.25 


340 ± 15 ± 70 


(+80, -50) 


0.75 


720 ± 25 ± 60 


(+100,-80) 


1.25 


700 ± 25 ± 50 


(+80, -70) 


1.75 


725 ± 25 ± 50 


(+90, -80) 




190 < W < 233 GeV 




-0.25 


305 ± 15 ± 60 


(+70, -50) 


0.25 


585 ± 20 ± 30 


(+70, -60) 


0.75 


555 ± 20 ± 50 


(+60, -40) 


1.25 


495 ± 20 ± 30 


(+50,-50) 


1.75 


500 ± 20 ± 40 


(+70, -50) 




233 < W < 277 GeV 




-0.75 


220 ± 15 ± 30 


(+60, -40) 


-0.25 


535 ± 20 ± 30 


(+80, -60) 


0.25 


500 ± 20 ± 40 


(+50, -40) 


0.75 


495 ± 20 ± 30 


(+50, -40) 


1.25 


475 ± 20 ± 30 


(+50, -40) 


1.75 


475 ± 20 ± 20 


(+50, -50) 



Table 3: Differential e + p cross section dajdif et for inclusive jet production integrated above 
Ej? 1 > 14 GeV in the kinematic region defined by Q 2 < 4 GeV 2 and in three regions of W for 
jets with a cone radius R = 1.0. Other details as in Table |l|. 
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dcr/drf et ± stat. ± syst. [pb] 


syst. i?jf*-scale [pb] 




134 < W < 190 GeV 




0.25 


230 ± 15 ± 30 


(+50, -40) 


0.75 


480 ± 20 ± 60 


(+60,-50) 


1.25 


450 ± 20 ± 30 


(+50,-50) 


1.75 


375 ± 15 ± 40 


(+50, -40) 




190 < W < 233 GeV 




-0.25 


215 ± 15 ± 40 


f + g0 -40) 


0.25 


400 ± 20 ± 10 


C+50 -40) 


0.75 


390 ± 15 ± 30 


C+40 -30) 


1.25 


320 ± 15 ± 30 


(+40, -30) 


1.75 


270 ± 15 ± 20 


(+30, -30) 




233 < W < 277 GeV 




-0.75 


165 ± 15 ± 30 


(+50, -30) 


-0.25 


385 ± 15 ± 30 


(+50, -40) 


0.25 


365 ± 15 ± 10 


(+40, -40) 


0.75 


320 ± 15 ± 40 


(+30,-20) 


1.25 


285 ± 15 ± 10 


(+30, -30) 


1.75 


255 ±15 ±10 


(+20, -30) 



Table 4: Differential e + p cross section dojdrf et for inclusive jet production integrated above 
Ej^ 1 > 14 GeV in the kinematic region defined by Q 2 < 4 GeV 2 and in three regions of W for 
jets with a cone radius R = 0.7. Other details as in Table IB. 



Cone radius 


a(R)± stat. ± syst. [pb] 


syst. -scale [pb] 


R = 0.5 


275 ± 10 ± 30 


(+40, -30) 


i? = 0.7 


385 ± 10 ± 20 


(+60,-50) 


1.0 


540 ± 15 ± 40 


(+90,-80) 



Table 5: e + p cross section <j(R) for inclusive jet production integrated above E^ 1 > 21 GeV 
and -0.5 < r] jet < 2 in the kinematic region defined by Q 2 < 4 GeV 2 and 134 < W < 277 
GeV. Other details as in Table |l|. 
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ZEUS 1994 




Figure 1: (a) The distribution of r tracks = E^\ racks j E 3 ^ l cal for the inclusive jet data sample with 
R = 1.0 (black dots) and as reproduced by the PYTHIA generator and detector simulation 
(histogram, normalised to the number of jets in the data); (b) the distribution of r di j et = 
E^l al (forward) /E^lai (central) for the dijet data sample (one jet in the forward region and the 
other in the central region) with R = 1.0 (black dots) and as reproduced by the PYTHIA 
generator and detector simulation (histogram, normalised to the number of jets in the data); 
(c) the values of the quantity (<rt racks >data / <r t racks>Mc) ~ 1 (circles) and {<r dijet > data / < 
fdijet>Mc) — 1 (squares). The shaded region displays the band of ±3% around zero. 
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Figure 2: Uncorrected transverse energy profiles for jets with > 14 GeV and R = 
1.0 and 0.7 as a function of the distance from the jet axis, Arj (integrated over \Af \ < tt/2), in 
three r]^ et regions (black dots). For comparison, PYTHIA and PYTHIA MI simulations includ- 
ing resolved plus direct processes are shown as the solid and dashed histograms, respectively. 
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Figure 3: Uncorrected transverse energy profiles for jets with > 21 GeV and R = 
1.0 and 0.7 as a function of the distance from the jet axis, Arj (integrated over \Af \ < tt/2), in 
three r]^ et regions (black dots). For comparison, PYTHIA and PYTHIA MI simulations includ- 
ing resolved plus direct processes are shown as the solid and dashed histograms, respectively. 
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Figure 4: Differential e + p cross section dajdrf^ for inclusive jet production integrated above 
Ej^ 1 from four different thresholds (E^ 1 > 14, 17, 21 and 25 GeV) in the kinematic region 
defined by Q 2 < 4 GeV 2 and 134 < W < 277 GeV for jets with a cone radius R = 1.0. 
The thick error bars represent the statistical errors of the data, and the thin error bars show 
the statistical errors and systematic uncertainties —not associated with the absolute energy 
scale of the jets— added in quadrature. The shaded bands display the uncertainty due to the 
absolute energy scale of the jets. For comparison, NLO calculations for three parametrisations 
of the photon parton distributions, = E^ 1 and for two different values of Rsep are shown: 
AFG R SE p = R (dot-dashed line), GRV-HO R SE p = R (dashed line), GS96 R SE p = R (thick 
solid line) and GS96 Rsep = 2i? (thin solid line). The values of Rsep used are indicated in 
parentheses. In all cases, the CTEQ4M proton parton distributions have been used. 
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Figure 5: Differential e + p cross section da/drf 
radius R = 0.7. Other details as in Figure f|. 
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Figure 6: Differential e + p cross section da/drf et for inclusive jet production integrated above 
Ej? 1 > 14 GeV in the kinematic region defined by Q 2 < 4 GeV 2 and in three regions of W: 134 
< W < 190 GeV (upper plot), 190 < W < 233 GeV (middle plot) and 233 < W < 277 GeV 
(lower plot) for jets with a cone radius R = 1.0. Other details as in Figure |j. 
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Figure 7: Differential e + p cross section da/drf** for inclusive jet production for jets with a cone 
radius R = 0.7. Other details as in Figure ||. 
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Figure 8: e + p cross section cr(R) as a function of the jet cone radius R for inclusive jet production 
integrated above E^ 1 > 21 GeV and —0.5 < if et < 2 in the kinematic region defined by 
Q 2 < 4 GeV 2 and 134 < W < 277 GeV. The thick error bars represent the statistical errors of 
the data, and the thin error bars show the statistical errors and systematic uncertainties —not 
associated with the absolute energy scale of the jets— added in quadrature. The shaded band 
displays the uncertainty due to the absolute energy scale of the jets. LO and NLO calculations 
using the GS96 (CTEQ4) parametrisations of the photon (proton) parton distributions and 
/j, = E^/A for two choices of the parameter Rsep are shown. The values of Rsep used are 
indicated in parentheses. The inset shows the calculations for a fixed value of Rsep — 2i2 and 
various choices of \i. 
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